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ABSTRACT 

The stability of nonaxisymmetric perturbations in differentially rotating astrophysical accretion disks 
is analyzed by fully incorporating the properties of shear flows. We verify the presence of discrete unstable 
eigenmodes with complex and pure imaginary eigenvalues, without any artificial disk edge boundaries, 
unlike Ogilvie & Pringle(1996)'s claim. By developing the mathematical theory of a non-self-adjoint 
system, we investigate the nonlocal behavior of eigenmodes in the vicinity of Alfven singularities at 
u>D = iWji, where ujd is the Doppler-shifted wave frequency and to a — k\\VA is the Alfven frequency. 
The structure of the spectrum of discrete eigenmodes is discussed and the magnetic field and wavenumber 
dependence of the growth rate are obtained. Exponentially growing modes are present even in a region 
where the local dispersion relation theory claims to have stable eigenvalues. The velocity field created by 
an eigenmode is obtained, which explains the anomalous angular momentum transport in the nonlinear 
stage of this stability. 

Subject headings: accretion, accretion disks — instabilities — MHD — plasmas 



1. INTRODUCTION 

Over the last several years, the presence of magnetic 
fields in a differentially rotating plasma has been pro- 
posed as a possible mechanism of accretion disk turbu- 
lence and its associated large anomalous angular momen- 
tum transport inside the disk. The presence of magnetic 
fields in a shear rotating gas cylinder ma kes the gas un- 



which depend strongly on the boundary conditions. 

In this paper, we concentrate our focus on the valid- 
ity of the linear analysis of nonaxisymmetric eigenmodes 
questioned by Ogilvie & Pringle(1996). The analysis is 
performed in the frame rotating with the local angular 
velo city, which is adopted in nonaxisymmetric mode ana l- 
ysis( Ogilvie fc Pringle 1996| ; Matsumoto fc Tajima 1995 ), 



stable against axisym metric perturbations ( Velikhov 1959 
Chandrasekhar 196l|). The normal mode analysis (Kumar. 



Coleman fc Kley 1994 ) of this local magnetoshearing in- 



stability showed the existence of unstable axisymmetric 
eigenmodes. The presence of this rob ust instability was 
re-recognized ( Balbus fc Hawley 1991) and confirmed by 
nonli near ideal MHD simulations flHawley fc Balbus 1991 



1992; Hawley, Gammie & Balbus 1995). The observational 
features of astrophysical accretion disks point to the need 
for large viscosity much beyond the one collisional mecha- 
nisms can yield. This rob ust instability has been invoked 
(Hawley & Balbus 1992) as a most pr omising candidate 
mech anism for the viscosity puzzle(e.g., Tajima & Shibata 
1997[ ). 

Numerical investigation of nonaxisymmetric magnetos- 
hearing modes has been carried out by adopt ing the 
shearing coordinates (e.g., Balbus fc Hawley 1992). Mat- 



sumoto & Tajima (1995) analyzed nonaxisymmetric non- 
local eigenmodes which are sandwiched by two Alfven sin- 
gularities around the corotational point and are not influ- 
enced by the disk edge boundaries and grow exponentially 
in time. These modes are distinct from the modes dis- 
cussed by Ogilvie & Pringle (1996), which are nonaxisym- 
metric modes contained within cylindrical boundaries and 



since eigenmodes evolve exponentially in time. The reso- 
lution of this question is important in theory of accretion 
disk Unless this magnetoshearing instability is a robust 
mode unaffected by the boundary conditions, the long 
search for candidate mechanisms of anomalous viscosity 
of accretion disks needs to be reopened. The criticism of 
Ogilvie & Pringle (1996) is interesting because it reflects 
the difficulty and the extreme mathematical and physical 
subtlety involved in the nature of this mode around the 
Alfven singularity. The problem is marred by the non- 
self-adjointness of the differential equation that describes 
eigenmodes, arising from the presence of shear flows. We 
know of no systematic mathematical theory on non-self- 
adjoint differential equations. Thus it takes us a develop- 
ment of mathematical and physical theory of such a system 
in order to understand the argument by Ogilvie & Pringle 
and to respond aptly to it. In the end such analysis has 
been developed and, assuringly, we find that the magne- 
toshearing instability is both robust and insensitive to the 
boundary condition, as we thought originally. 

In §2, we derive the wave equation in a differentially 
rotating magnetized disk based on the analysis of Mat- 
sumoto & Tajima (1995), and basic properties of the wave 
equation are discussed. We show that discrete nonaxisym- 
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metric eigenmodes exist, which are buffeted by the pair of 
Alfven singularities where the Dopplcr-shifted wave fre- 
quency equals the Alfven frequency. Our analysis finds 
that the eigenmodes oscillate indefinitely in the vicinity 
of the Alfven singular points when the eigenvalue is real, 
whereas the eigenmodes are regular when the eigenvalue 
is complex. Numerical calculation of eigenmodes with the 
Alfven frequency and wavenumber dependence is discussed 
in §3. We compare the results with the local dispersion 
relation, and show these eigenmodes are discrete. Astro- 
physical implications and conclusions are discussed in §4. 

2. ANALYTICAL PROPERTIES OF NON-SELF-ADJOINT 
EQUATION NEAR ALFVEN SINGULARITY 

We consider the MHD stability of magnetoshearing 
modes in the co-rotating frame of the fluid. The basic 
ideal MHD equations in the frame rotating with angular 
velocity f2 are 



where u>d is the Doppler-shifted frequency, 



lud = OJ + - flkyX, 
and to a is the Alfven frequency, 



(7) 



(8) 



The initial condition enters through the source function 
T(x,iv). 

The wave equation is derived by expressing the homo- 
geneous part of equation (^|) in terms of the normalized 
radial coordinate 

£ = ¥^ (9) 
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+ 2v x n + (Q, x r) x fi, (1) 



— =V x (v x B), 
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where g is the gravitational acceleration and r is the po- 
sition vector. We assume incompressibility for simplicity, 



V • v = 0. 



(3) 



We also ignore self gravity, which is not essential for the 
magnetoshearing instability. 

We use the local Cartesian coordinates (x, y, z) in the ro- 
tating frame where the ir-axis is in the radial direction, the 
y-axis in the azimuthal direction, and the z-axis parallel to 
f2. The uniform velocity shear v y = — (3Q/2)x is assumed 
for the Keplerian disk, where x = is the local co-rotating 
radial position. The wave equation is derived by lineariz- 
ing the basic equations around the equilibrium state and 
assuming solution of the form 4>(x, i)exp[i(fc a y + k z z)]. In 
the unperturbed state, the density, pressure and mag- 
netic field are assumed to be uniform. The assumption 
of v x = v z = B x — in the unperturbed state yields the 
unperturbed momentum equation, 



g + 2v Q x Jl + (SI x r) x ft = 0. 



(4) 



Next, the Laplace transform of the perturbation, 
(x,li), is employed, 



(j>(x,u) = / dt(j>e h 
Jo 



(5) 



Substitution of the Laplace transformed momentum and 
indu ction equations into the c ontinuity equation yields 
(see Matsumoto fc Tajima 1995 for detail) the initial value 
equation 
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v x = T(x,co), (6) 
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where the ratio of the squares of the azimuthal and vertical 
wavenumber is defined as 



At 



(11) 



Unstable eigenmodes may exist when the solution satis- 
fies the boundary condition 



lim v T 







(12) 



in the upper half of the complex cj-plane reference. This 
boundary condition makes our eigenmodes distinct from 
the modes found by Ogilvie & Pringle (1996), which are 
confined in rigid cylindrical boundaries and strongly de- 
pendent on the boundary condition. In order to have 
an overall angular momentum transport across the entire 
disk, it is imperative to have unstable modes within the 
disk, not just on the boundaries of the disk. To investigate 
the interior of the accretion disk, eigenmodes should not 
depend on the edge boundary conditions. The eigenmodes 
which arise from finite boundaries may contribute to the 
angular momentum only near the edge of the disk. Our 
mode, which is not affected by the edge, grows whenever 
the eigenfunction is located between two Doppler-shifted 
Alfven points. Since the positions of those points are de- 
termined in co-rotating frame, the origin of the frame can 
be anywhere in the disk. The boundary condition then 
assures that the momentum transport resulted from the 
superposition of the growing eigenmodes which can occur 
throughout the accretion disk. 

The proper boundary condition interior of the disk can 
be easily examined by inspecting the aysmptotic form of 
the basic differential equation (0) of the system. This in- 
dicates that the leading radial dependence of eq. (||) leads 
to the exponential decay away from the co-rotation point 
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and toward the Alfven singularity. This mathematics is, 
of course, most reasonable and physical as well, because 
the i nstability energy is provided within the two Alfven 



singular layers to the mode, which dissipates the energy 



at (or near) the singularities. 

Since the wave equation ( [To| ) is not self-adjoint due to 
the existence of the flow shear, the square of the eigenvalue 
lo 2 is not guaranteed to be real, which any self-adjoint sys- 
tem always satisfies. The fact that the eigenvalue is not 
pure real or imaginary but in general complex prevents us 
from applying the Strum-Liouville theory to this system. 
Note that if k y = 0, equation ffl) is reduced to self-adjoint, 



dx 2 
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(13) 



and its eigenvalue constitutes the Alfven continuum. 
The analysis of this mode has been done by Chan- 
drasekhar(1961). We now assume k y ^ 0. 

To the best of our knowledge, no theory for non-self- 
adjoint operators exists, and we have to investigate the 
properties of non-self-adjoint systems in general. We find 
that the eigenmodes possess certain symmetry properties, 
which are originated in the character of the differential op- 
erator and radial symmetry. Since the differential operator 
D(lo, £) is invariant under the operation (w, £) — » (— lo, — £) 
because of radial symmetry, the eigenfunction v x (lo, £) is 
also the eigenfunction of D{— lo, — £), 



D(-w,-0«x(w,0 = 0. 



(14) 



Note that the operation £ — > — £ is the same as changing 
the direction of the rotation f2 — > —51. Taking conjugate 
of the equation (Oh and denoting — £ to £ yields 



D(-lo*,$v x (lo,-0 
= U(-w*,0«x(-w*,0 = > 



(15) 



where <p* means the complex conjugate of (f>. Thus, if lo is 
an unstable eigenvalue of the equation (ppj), — lo* is another 
unstable eigenvalue whose eigenfunction v x (~lo* ,£) satis- 
fies the relation fllq ) . Especially, when lo is pure imaginary, 
i.e., — lo* = lo, the real part of eigenfunction is symmet- 
ric and the imaginary part antisymmetric with respect to 

e = o, 

«i(w,-0 = *»(w,0- (16) 

These properties of the non-self-adjoint operator indicate 
that this system has complex eigenvalues in general, which 
differs from a self-adjoint system, and one unstable eigen- 
value has another unstable and two stable companions. 
This symmetry of non-self-adjoint system and the compar- 
ison between self-adjoint and non-self-adjoint eigenvalues 
in complex-w plane is shown in Figure 1. 

Next, we look for the solution of the wave equation ([Tof). 
The boundary conditions for these ideal MHD modes we 
are interested in are that the eigenmode resides and is dif- 
ferentiable around the corotational point and sandwiched 
by a pair of the Alfven singularities, and decays toward 
£ — > ±oo. Note that these physical boundary condi- 
tions preclude t he usual Kelvin- Hc l mholtz(K-H) instabi l- 
ity eigenmodes (Tajima et al. 1991; Noguchi et al. 1998| ) . 
The boundary conditions generic to the K-H mode is to 
match the dissipative structure at the corotational point. 



The properties of these eigenmodes may be investigated 
by analyzing the solution around the spatial co ordinates 
of interest by using the Frobenius expansion (Morse & 
Feshbach, 1953| ), in particular around the Alfven singular 



points. It is of particular significance to examine analytical 
properties of the Frobenius expression around the Alfven 
singularities in order to examine the question leveled by 
Ogilvie & Pringle (1996) against the eigenmodes obtained 
by Matsumoto & Tajima(1995). An appropriate, compat- 
ible numerical method for these eigenmodes is, therefore, 
the shooting method that starts from an exponentially de- 
caying functional form at £ = ±oo and shoots toward the 
corotating point where two sides of the function should 
smoothly (differentiably) match. (This will be closely ex- 
amined in the next section.) 

We expand v x in equation (ll^) using, 



n=0 



M-U±) 



n+s 



in the vicinity of the Alfven singularities lojj 
£ = £a±, where £ A ± is defined by 

LO 

U± = ±1 

LO A 

= ±1 — lo' — itoL 



(17) 



±loa or 



(18) 



and Lo'rico'i) is the real(imaginary) part of lo/loa- Applying 
the Frobenius method in the vicinity of the regular singu- 
lar points £_= £,a±, the indicial equation for the exponent 
s for Eq. (AG) is given by 



s z + 



2 T - uj[ 
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4 T 2iw^ 
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0. (19) 



Note that the eigenfunction is irregular and oscillates in- 
definitely at £ = £a± whenever s has imaginary compo- 
nent, provided Re(s) is nonpositive integer, 



«x = (£ - U±r 4s) exp [Im(s) log |£ - ^±|] 



(20) 



The singular points are on the real axis if and only if the 
eigenvalue lo is real. In general, the indices are given by 



u'M ±* 3i) 



2(2 T *i^)(l±* H) 



±4 1- 



16 (4T*2i^-^')(lT*i^-) 5 



9q 



<04±* 3i) 2 



(21) 



where the sign ±* indicates that we take upper sign when 
£ = ^a+ , and lower sign when £ = £a- j respectively. Since 
s is complex, the solution is not analytic at the Alfven sin- 
gular points. 

Now, we investigate some special cases. First, when the 
eigenvalue is pure real, i.e., lo[ — 0, the indices are purely 
imaginary, 

S = ± ^ W 

the eigenfunction rapidly oscillates and is indefinite. Note 
that the boundary condition in this case becomes special 
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Fig. 1. — Symmetry of a nonself-adjoint system and the comparison between self-adjoint and nonself-adjoint eigenvalues in complex-o; 
plane. Discrete modes are denoted by crosses. The eigenvalue of the self-adjoint operator should be purely real or imaginary, and the eigen- 
functions are predictable by Strum-Liouville theory. However, there is no such restriction for the eigenvalue of the nonself-adjoint operator, 
and u),u)* , — lu and — ui* make a group of solutions. 



in that we need to shoot outward from £ = £ c . We can no 
longer shoot from the outer to the inner region. As the 
eigenmodes oscillate indefinitely around the Alfven singu- 
larity, we find that the eigenvalue becomes continuous in 
this pure real case to form the AlfVen continuum. Second, 
when oj[ <C 1, which is similar to the first case but now 
the singular points are not on the real axis, the solutions 
are 

" 2 . . 3 



s = ± 



3^32^. 



(23) 



whose solution oscillates limited times in the vicinity of the 
singular points on the real axis, but not indefinitely, and 
the solution is regular on the real axis. Third, when the 
perturbation is nearly toroidal (k y k z ,q^> 1), as in the 
accretion disks far from their source, and the eigenvalue is 
complex, s is given by 



u'M ±* 3i) 



2(2 T *i^)(lT* H) 



± 1 



8 (4T*2i^-a,f)(lT*M) s 



9q 



<H ±* 3i) 2 



(24) 



When the perturbation is pure toroidal, the second term 
in the square brackets in equation ( pif ) vanishes and s is 
given by 

4(4±*3i) 
' (2T*i^)(lT*i-0- { ' 

The eigenfunction corresponding to s — is regular, 
and the latter one diverges logarithmically at the singu- 
lar point, when ui is real. Finally, when the perturbation 
is nearly poloidal^ <fc z ,g< 1), as may occur in accre- 
tion disks close to their source, we find s as 



± 



2i ^4t'2^-^ 
3^q 2 T* H 



(26) 



which reduces to the roots of the first case when the eigen- 
value is pure real. 

The exponent s at the corotation point lud — (or 
£ c = lo/oja) is given by 



1 



2K 2 - 1) I 



3 ± \ u>. 
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(27) 



When uj'^ > 7 + 4\/2, the eigenvalue is regular at the 
corotation point, since both indices are real and positive. 
When < c^ 2 < 1, indices are still real, and one of them 
is positive. This is the solution that is consistent with the 
matching condition at £ = of the shooting method dis- 
cussed in §3. When 1 < cj- 2 < 7 + 4v / 2, the corotation 
point is singluar since two indices are real and negative in 
the region 1 < lj' 2 < 7 — 4\/2, and complex in the region 
7-4V2 < wf < 7+4V2. When wf = 1, the eigenfunction 
has irregular singularity at the corotation point. Again, 
even though the eigenfunction is irregular at the corota- 
tion point, the physical eigenmodes on the real £-axis is 
regular. Moreover, since all the coefficients of differential 
equation (10) are real at the corotation point even when 
the eigenvalue is complex, the eigenfunction should be real 
at f = £ c . 

Although irregular in the vicinity of the Alfven singular- 
ities or the corotation point in most cases, eigenfunctions 
are not irregular in the physical sense unless the singular- 
ities are on the real axis. Instead, the oscillatory behav- 
ior and amplitude of the eigenfunction around the Alfven 
singularities directly reflects the physical eigenfunction be- 
havior on the real £-axis, especially if is small, i.e., mag- 
netic field is strong. 

When the eigenvalue oj and the index s are both com- 
plex, the eigenfunction oscillates indefinitely in the vicinity 
of the singularit y d ue to the imaginary component of the 
index s(see eq. |]20| ). In this case, the physical eigenfunc- 
tion on the real £-axis also oscillates very rapidly in the 
vicinity of the point which is the projection of the complex 
singular point to the real £-axis, but the physical eigen- 
function oscillates only finite times because the projected 
point is not a singular point. 

When the real component of s is negative, the eigenfunc- 
tion diverges at the singularity(eq. Q ) ■ The amplitude of 
the physical eigenfunction on the real £-axis is large at the 
projected singular point on the real £-axis. However, since 
the projected point is not a singular point, the eigenfunc- 
tion does not diverge at this point. It is also clear that the 
eigenvalue is regular even if the singular point is a branch 
point, since we can choose the branch cut of the eigenvalue 
without crossing the real £-axis. 

If the eigenvalue is real, the singularities are on the real 
£-axis and eigenfunction is irregular in the physical sense. 
We will discuss pure real eigenvalue cases in §3. 
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3. ROBUSTLY UNSTABLE MAGNETOSHEARING 
EIGENMODES 

The eigenvalues of the wave equation ( pX| ) are calculated 
numerically by the shooting method with the boundary 
condition discussed in §2. Since the equation ( |To| ) may 
have three singularities (corotation point and Alfven reso- 
nances), we choose complex initial value to avoid Alfven 
singularities on real axis, and integrate(" shoot") the equa- 
tion (Il0|) on the real £-axis from the left and right asymp- 
totic boundaries (which are for removed from the any 
of particular singular behavior) to the corotation point, 
where their value and first derivative of the eigenmode are 
to be matched. If they are not matched, we change the 
eigenvalue appropriately until we can match. This itera- 
tive method is generally called the " shooting method" for 
eigenvalue problems. Spatial steps of the integration are 
smaller in the vicinity of the point where the Alfven point 
is projected on the real £-axis than other regions, in case 
the eigenvalue is almost real but still complex. By using 
the Newton method to decrease errors of a trial function, 
we obtain higher accuracy and faster convergence than the 
previous shooting codes. This allows us to search for sub- 
tle singular and regular eigenfunctions over a wide range 
of parameter values. 

In order to satisfy the boundary condition (|l2]), we im- 
pose v' x /v x = k± at the numerical boundaries £ = ±10 
(corresponding to the artificial infinity) , where k± are the 
negative and positive solutions of the quadratic equation 
given by inserting the functional form v x = exp(k±x) into 
the equation (px|). On the numerical boundaries, the lead- 
ing term of the equation (ra) is the first term of the coef- 
ficient of v x , which is of the order of 1, and other terms 
are of the order of 0.01 or lower in our calculation. Our 
assumption for the boundary condition is justified as far 
as these estimations are valid. 

Figure 2 shows examples of eigenfunctions v x obtained 
by our shooting code when lua — 0.1ft and q — 0.01. The 
solid and dashed curves represent the real and imaginary 
parts of the eigenfunction respectively. Figure 2a is for the 
fundamental pure imaginary eigenvalue and Fig. 2b is for 
the complex eigenvalue. Since the eigenvalue of Fig. 2a is 
pure imaginary, the real part of the eigenfunction is sym- 
metric and the imaginary part antisymmetric with respect 
to £ = 0, which is consistent with the equation (|l6|). Figure 
2b is the eigenfunction with a complex eigenvalue, which 
makes a pair with the eigenvalue — lu* whose eigenfunction 
is derived from the relation (|l5|). These eigenfunctions 
are confined between two Alfven singularities located at 
£ = £a± ~ ±1, and they are real at £ = £ c . 

Figure 3 shows the distribution of eigenvalues in the up- 
per complex w-plane when lua — 0.01 and q = 0.01. It 
shows only the eigenvalues in the region Re(w) > 0, and 
all the complex eigenvalues have a paired eigenvalue — lu* 
in the region Re(w) < 0. It is obvious that this non-self- 
adjoint system has complex eigenvalues, which never ap- 
pear in a self-adjoint system (see Fig. 1). There are only 
two pure imaginary eigenvalues, which will be shown to 
merge by changing lua and q. We find that complex eigen- 
values, which Matsumoto & Tajima(1995) did not find, 
exist and have smaller imaginary part and grow slower in 
time than the fundamental eigenmode. 

Figure 4 shows the dependence of unstable eigenvalues 



on lua when q = 0.01. The solid (dashed) curves show 
the imaginary(real) part of the eigenvalues. When lua is 
small, there exist two purely growing modes, which merge 
at lua ~ 0.66ft and form complex eigenvalues. These 
modes were found in Matsumoto & Tajima (1995) and 
the qualitative properties of this mode are about the same 
as found in Matsumoto & Tajima. However, the merging 
point is slightly greater than the earlier value and, more 
significantly, the growth rate does not decay significantly, 
even beyond lua = 1.5840, where it was calculated to van- 
ish in Matsumoto et al.. The fundamental mode acquires 
its maximum growth rate just before it merges with an- 
other (secondary) pure imaginary mode. In addition, we 
find two new complex modes. These complex eigenmodes 
are also shown in Fig. 4, and the growth rate of all these 
modes saturate to lu ~ 0.15fti, the same saturation value 
for the fundamental mode with increasing lua- 

Figure 5 shows the dependence of eigenvalues on q when 
lua = 0.01ft (Fig. 5a) and 0.66ft (Fig. 5b). The funda- 
mental^) and secondary(S) pure imaginary modes and 
two complex modes(l, 2) are shown in Fig. 5, which cor- 
respond to the eigenmodes in Fig. 4, when q = 0.01. 
Eigenvalues calculated from local analysis are also shown 
in Fig. 5, which we discuss later. Two pure imaginary 
modes are always distinct when lua = 0.01ft (two upper 
modes in Fig. 5a). However, when lua = 0.66ft, these two 
modes merge and become complex at log 10 q ~ — 1.8 and 
split to become pure imaginary again at log 10 q ~ 1.4. The 
eigenvalues of the other two complex modes become pure 
imaginary when q exceeds a certain value (logg = —1.25 
for lua = 0.01ft, logq = -0.75 for lu a = 0.66ft), and the 
growth rate for those modes saturates with increasing q. 

Next, we compare the nonlocal eigenfunction results 
with the local(Fourier) dispersion relation. By replac- 
ing d/dx in equation (^) with a constant \k x around 
x = and assuming that the unperturbed magnetic field 
is toroida^-B^ = B z =0), the local solution in the regime 



M ~ LUA « OUJ = y/Wkl/ikl + kj+ kl) 

is (Matsumoto & Tajima 1995) 
3 



1 



2 q± 2 



(q-2)(q--) 



2 

LU A - 



(28) 



This local dispersion relation ([28|) shows that pure real 
eigenmodes appear in the region q < |, pure imaginary 
eigenmodes in q > 2, and complex in | < q < 2. How- 
ever, when q is small, i.e., the perturbation is almost par- 
allel to the magnetic field, nonlocal eigenmodes are un- 
stable in both lua = 0.01ft and 0.66ft (see Figs. 5a, b), 
and the growth rate of each mode does not have strong 
dependence on q. We conclude that replacing d/dx by 
a single wavenumber k x is invalid for these modes since 
such eigenmodes oscillate very rapidly in the vicinity of 
the Alfven points in a pronounced fashion(see Fig. 2b). 
In other words, the spatial variation of the wavenumber in 
the radial direction is essential for the modal analysis of the 
magnetoshear instability. Radial dependency of the wave- 
length also prevents us from applying the WKB method 
to this model. The WKB method requires the wave- 
length of the eigenmodes L e is smaller than the shear scale 
length L s (L e /L s <C 1), which may be satisfied around 
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Fig. 2. — Examples of the eigenfunctions of the nonaxisymmetric magnetoshearing instability in a Keplarian disk. Solid curve, real part 
of eigenfunction; dashed curve, imaginary part of eigenfunction. The model parameters are u)a = 0.1C2 and q = 0.01. (a) Fundamental pure 
imaginary mode (eigenvalue is u = 0.00785f2i);(£>) Complex mode(eigenvalue is u) = (0. 00109 + 0. 00039i)f2). The eigenfunctions are sandwiched 
by two Alfven singularities £ = ±1 — uj/wa- 
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Fig. 3. — Distribution of unstable eigenvalues of the magnetoshearing instability in the upper complex w-plane when lua = 0.1S1 and 
q = 0.01. Eigenvalues with Re(o;) < are shown, and all the complex eigenvalues have a paired unstable eigenvalue — lu* in the region 
Re(oj) < 0. There exist only two pure imaginary eigenmodes and many complex eigenmodes, which are not allowed to exist in self-adjoint 
system. 
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Fig. 4. — The Alfven frequency (u>a = k\\VA) dependence of eigenvalues of magnetoshearing instability when q = 0.01. The dashed curves 
and solid curves show the growth rate lm(uj) and the real frequency Re(u;), respectively. The fundamental(F) and secondary(S) pure imaginary 
eigenmodes and two complex modes(l, 2) are shown, which correspond to the eigenmodes labeled F, S, 1 and 2 in Fig. 3, respectively, when 
uia/Q = 0.01. Two pure imaginary eigenvalues merge at loa ~ 0.66O to form complex eigenvalues. The growth rate of all four eigenvalues 
saturates to u> ~ 0.15f2i with increasing lua- 




the corotational and Alfven singularities, but the wave- 
length is comparable to the shear scale length in other 
regions(L e /L s ~ 1). 



To show that these eigenmodes are discrete, let us first 
show the existence of the Alfven continuum on the real 
w-axis. The wave equation (10) has a solution for any 
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Fig. 5. — The q(=ky/k^) dependence of eigenvalues of magnetoshearing instability. The dashed curves and solid curves show the growth 
rate Im(oj) and Re (u>), respectively. The fundamental(F) and secondary(S) pure imaginary modes and two comlex modes(l, 2) are shown, 
which correspond to the eigenmodes labeled F, S, 1 and 2 in Fig. 4 when q = 0.01. Eigenvalues calculated from local mode analysis are 
also shown, (a) u>a = 0.01S1. Two pure imaginary eigenmodes(F, S) are always distinct, and the complex eigenmodes(l, 2) become pure 
imaginary, whose growth rate saturates, with increasing q. (b) wa = 0.66C Two pure imaginary eigenvalues(F, S) merge at log 10 q ~ —1.8 
to form complex eigenvalue, and split again to become imaginary at log 10 q ~ —1.4. The complex eigenmodes(l, 2) becomes pure imaginary 
with increasing q, and the growth rate saturates. In both cases, local modes are stable even in a region where nonlocal modes are unstable. 



real u> for which ujfj = lo\ for some x. It follows that 
the spectrum of this mode is continuous, and the Alfven 
continuum extends to the entire real u> by choosing some 
-B, k y and k z , which is different from the model chosen by 
Ogilvie & Pringle (1996) in which the Alfven continuum 
is restricted by the boundary condition. The eigenmodes 
with the pure imaginary eigenvalue that we have shown 
are obviously not in this class. 

When B — 0, the Kelvin- Hclmholtz modes can be de- 
rived, which are stable in accretion disks. In this limit, 
equation (0) reduces to 



d 2 v x 
dx 2 



-Hi 



2 1,21 



v x = 0, 



(29) 



and when k z = 0, it has a simple solution v x = 
exp[—k y \x\], which satisfies the boundary condition (|l2|). 
The first derivative of this class of solutions is discontin- 
uous at x = 0, which vanishes with introducing dissipa- 
tion. The Kelvin-Helmholtz instability also has continu- 
ous eigenvalues, but this class of solutions is eliminated in 
our calculation because of the matching condition of the 
shooting method, which requires the eigenfunction and its 
first derivative to be continuous. We search for eigenvalues 
by choosing initial values in the region < Lu r /ujA < 1 and 
< uJi/uJA < 1- We find that of the initial values converge 
to one of the eigenvalues in Fig. 3 and we conclude that 
all of the eigenmodes are discrete. 

In Fig. 6 we show another eigenmode whose eigen- 
value gradually becomes real with increasing q, when 
oja = 0.6611. However, when the eigenvalue is real, we 
have already shown that the index of the eigenfunction s 
is pure imaginary in the vicinity of the Alfven singularities 
(Eq. (22)) and that the eigenfunction has the form 



= exp [i|s|log(£-£4±)] 



(30) 



from equation (p0|). Such eigenfunctions oscillate indef- 
initely toward the singularity and the function can take 
any value between — 1 < v x < 1. This indicates that the 
function in the inner region £a_ < £ < £a+ and outer 



region £ < ^a-,S, < £a+ is discontinuous at the Alfven 
singularities. Thus the boundary condition for the con- 
tinua cannot be that of shooting from the outside |£| = oo 
toward the inside, but we should shoot from inside toward 
the singularities. 

Figure 7 shows an example of eigenfunction in the inner 
region when oja = O.Olfi, w = 0.002S1 and q — 0.01 (see 
eq. (22) and arguments for detail of this mode). How- 
ever, the eigenmode in Fig. 6 is continuous even at the 
Alfven singular points, since the integration by a finite 
spatial step brings in an effective dissipation, which is not 
the case for the pure real eigenvalue. Instead, the eigen- 
mode becomes continuous because of the numerical dissi- 
pation. Although this numerical eigenfunction is different 
from the theoretical eigenfunction beyond the passage of 
the singularity, the fact of continua remains the same for 
two different reasons. It should be pointed out that in real 
physical situation there always exists a dissipation even for 
a nearly ideal MHD system. The dissipation prevents the 
eigenmode from blowing up on the Alfven singular points, 
keeping the energy of the eigenmode finite. The numeri- 
cal dissipation affects this eigenmode in the same manner 
mathematically as the real dissipation does, by passing 
oscillations through the singularity barrier and avaraging 
oscillation in a finite spatial step. Thus the numerically 
obtained eigenmode, though different from theoretically 
expected continua, may be regarded as realistic. 

Finally, we describe the physical behavior of the eigen- 
modes in accretion disks. The expression of v y in terms of 
v x is derived from the continuity equation^, 



q d 
k~ z dt 



1 



+ q 



2K - wi) 



( 3 4- 


-)] 







(31) 



Figure 8a shows v y calculated from the fundamental pure 
imaginary mode (Fig. 2a) and the velocity field created in 
the xy-pl&ne by the fundamental eigenmode is shown in 
Fig. 8b. The eigenfunction of v y is also trapped between 
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Fig. 6. — The example of a complex eigenmode that becomes pure real value with increasing q when u)a = 0.66Q. The dashed curve and 
solid curve show the growth rate Im(o;) and the real frequency Re(aj), respectively. The eigenvalue is complex when q is small, which becomes 
real with increasing lua. 
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Fig. 7. — The mode with a pure real frequency that has the Alfven singularities at ui = ±U/i, where the mode energy evaneces. 
eigenfunction oscillates indefinitely towards ui = uij\, and also towards uj = —uja with small amplitude. 
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Fig. 8. — An example of v y calculated from the fundamental pure imaginary mode (Fig. 2a) and the velocity field in xy plane by the 
fundamental mode. The solid curve and dashed curves show the real and imaginary part of v v , respectively, (a) v y is almost out of phase 
with respect to v x . (b) Velocity field is created in xy-plane by the fundamental mode. Vortices are created between two Alfven singularities, 
and they will overlap with other eigenmodes excited at various ^-positions to expand the unstable region. 



two Alfven singularities. Since both v x and v y are almost 
out of phase with each other, the velocity field created by 
the fundamental eigenmode consists of vorti ces in xy-plaiie 
which arc the seeds of nonlinear instability (Matsumoto & 
Tajima 1995). Note that since we use the frame rotating 
with angular velocity f2, there is no unique origin x in the 
xy-plane. Thus such unstable eigenmodes excited at vari- 



ous x-positions will overlap with each other to expand the 
unstable region in the x-direction. 

4. SUMMARY 

We have shown the existence of the discrete unstable 
nonaxisymmetric magnetoshearing instability eigenmodes. 
Since we assumed the exponentially decaying boundary 
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condition (Eq. |12] ) for the radial component of velocity, 
our modes are independent of the effect of the boundary 
condition. In our accretion disk theory, this robust in- 
stability occurs without any unrealistic disk edge bound- 
ary condition in infinite linear shear flow. The scale- 
length of a single eigenmode in the x direction Ax is de- 
termined by the local strength of magnetic field, the di- 
rection and amplitude of the wavenumber and the mag- 
nitude of the angular velocity (see Eqs. @ and J8|), 
since the eigenfunction is buffeted by the Alfven singular 
points up = ±uja- If the magnetic field is pure toroidal, 
Ax = 2w A /3ttk y = 2v A /3n ~ 2/3(v A /C s )H, where C s 
is the sound speed and H = C s /fi is the thickness of the 
disk. When v A *C C s , the mode is localized in the radial 
direction with the scalelength smaller than the thickness of 
the disk. If the magnetic fie ld has azimuthal component , 
Ax is propotional to fc^/fc^ ( [Matsumoto fc Tajima 1995j ). 
In both cases, our infinite boundary condition is sufficient 
if the scalelength of the eigenmode is smaller than H . The 
curvature of the magnetic field is also small if Ax <C H . 
However, for nearly axisymmetric perturbations (fey <C fen), 
the eigenmode have a large radial scalelength, and the in- 
finite boundary condition may not valid. Density gradient 
and geometrical effects become important in this case. 

Our result of the growth rate of unstable modes agrees 
with that of Matsumoto & Tajima (1995) in the region 
oj a < f2. We have found complex eigenvalues with smaller 
growth rates than the fundamental pure imaginary eigen- 
mode. When uj a is larger than fi, two pure imaginary 
eigenmodes merge, the results of which is the same as 
found in Matsumoto & Tajima. However, our result shows 
that the growth rate saturates with increasing ui A . This 
indicates that the accretion disk is unstable even if the 
Alfven frequency is comparable to the angular velocity, a 
case of strong magnetic fields. 

The comparison of the nonlocal and local dispersion 
relations demonstrates where and how the local Fourier 
mode approximation fails to be accurate for this nonax- 
isymmetric instability. The wavenumber dependence of 
the eigenvalue shows that the nonlocal modes are unstable 
even in the region g = ^/^ < 1, where the local disper- 
sion relation has only a stable solution. Furthermore, we 
have found that two pure imaginary eigenvalues merge to 
be complex and split into two pure imaginary again with 
increasing q in a region where the solutions of the local 
dispersion relation are purely real. Overall, as Fig. 5 indi- 
cates, the discrepancy of the local theory from the correct 
nonlocal theory amounts to not just a quantitative level 
but a qualitative deviation. 

We have developed the mathematical and physical the- 
ory of a system of nonself-adjoint differential equations for 
the first time. In astrophysics, nonself-adjointness always 



appears whenever there is a shear flow, which is mathemat- 
ically unsolved so far. We find the relationship of complex 
eigenvalues, which a self-adjoint system does not have, by 
a general approach to the nonself-adjoint system. A pair of 
eigenvalues uj and —lu* relate to each other in our model, 
since our model is symmetric with respect to the x-axis. 
Even if there is no such symmetry, we conclude that four 
eigenmodes make up a group in the nonself-adjoint system 
in general. 

Although our model of the nonaxisymmetric mode has 
been linear in Cartesian coordinates and ignored the effect 
of diffusion, analysis in Section 3 suggests how the eigen- 
mode grows to enter a nonlinear stage, and how it explains 
momentum transport in accretion disks. Since there is no 
particular sense in the radial direction in Cartesian coor- 
dinates, there exists no specific calibration of momentum 
transport in the linear stage. The eigenmodes can be ex- 
cited, however, between any pair of Alfven singularities 
in the radial direction (the x-axis) and create vortices in 
the disk plane (the xy-plane), as shown in Fig. 8b. The 
eigenmode with the fundamental eigenvalue dominates in 
time for a given radial co-rotation point. For another 
(arbitrary) co-rotation point, the same applies. These 
eigenmodes can overlap with each other to form greater 
vortices. In this stage, the nonlinear effect gives rise to 
anomalous magnetic viscosity that underlies the momen- 
tum transport needed to explain astrophysical disks. Mat- 
sumoto & Tajima demosntrated non-linear evolution of the 
eigenmode by three-dimentional MHD simulation with the 
shearing-box model. The overlap of eigenmodes excited at 
various xs' was shown. They also calculated the magnetic 
viscosity parameter 

(6B x 6B y ) (SB") „, ^ 1 

where the notation (5B x 5B y ) etc., denotes the spatial 
avarage. They found that when the poloidal field is domi- 
nant, the magnetic viscosity is ~ O(0. 1), which corre- 
sponds to a in dwarf novae during the bursting phase. 

We conclude that the results of Matsumoto & Tajima 
are correct and that the robust mode in magnetized ac- 
cretion disks is of the magnetoshearing origin. This mode 
should be dominant in nonlinear theory and our linear 
analysis supports the results from the three-dimensional 
simulation in Matsumoto & Tajima, which explained 
anomalous momentum transport in accretion disks. 

The work was supported by the US Department of En- 
ergy and NSF ATM 98-15809. TT is currently on leave at 
LLNL. 



APPENDIX 



EXISITENCE OF THE LOCALIZED GROWING MODE IN THE INFINITE SHEARING FLOW 

Ogilvie & Pringle's argument of the non-existence of localized growing modes ( pgilvie fc Pringlc 1996 , Appendix C), is 
incorrect for the following reasons. 

First, we rewrite the integrated wave equation ( Ogilvie fc Pringle 1996| , [C8]) without introducing a parameter A for 
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simplicity, 



where k 2 



+ OC 



dy_ 

dx 



(k 2 -f)\yf 



dx = 0, 



(Al) 



k 2 

Ky 



k 2 and the eigenfunction y(x) is assumed to decay exponentially as |xj — > oo, as we assume for the 



boundary condition. They also assume a uniform magnetic field, which makes equation (Al) qualitatively equivalent to 
our wave equation (|Io|). The function f(x) is given by 



/(*) = 



(kyU 



[(iu>i — kyu'x) 2 
4f2 2 (iwj — kyu'x) 



,212 



/l, „ 



2Q,u' 



[{iuji - kyu'x) 2 - uj\] 



[(iuii — kyu'x) 2 — uj 2 



(A2) 



where u(x ) is the shearing velocity, and f(x) has two Alfven singular points in the lower half-plane. We s how here that 



equation (Al) can be satisfied because of those two Alfven singular points, which makes the integral (Al) equal to zero 
if the complex frequency is properly chosen. The positions of the Alfven singular points are given by 



\0J l ± LP A 



(A3) 



where is the imaginary part of the eigenfrequency, and k y is the wavenumber parallel to the shear flow, respectively. 
Note that the total integrand of the equation (Al) becomes positive and decays to zero as \x\ — > oo since f(x) decays to 
zero as |x| — > oo. 

Now, if k 2 > \f(x)\ for any x, the integrand is always positive and equation (Al) can not be satisfied with any 
eigenfrequency. However, if the singularity is near enough to or on the real axis, becomes large enough to satisfy 

\f(x)\ > k 2 in the vicinity of the singularities and the total integral can take any value. Note that the real part of the 
eigenfrequency does not affect the position of the singularity. Furthermore, if the eigenvalue is complex, as we calculated 
in §3, the integral can equal to zero without crossing any branch cut of the singularities. The refore, localized growing 
modes should exist by choosing eig enfr eq uencies which makes the integral of the equation (Al) equal to zero. In other 
words, the contour of the integral (Al) is not free to pass around the circle in the upper half plane, since the complex 
conjugate of the eigenfunction y* has two singularities in the upper half plane with branch cuts, which any half circle in 
the upper half plane with sufficient large radius should cross. 

Thus the statement "The integral relation (C8) therefore cannnot be satisfied" in line 25 of p. 164 of Ogilvie & Pringle 
was incorrect, and our calculation is valid. 
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